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Skyrmions are topologically protected entities in magnetic materials which have the potential to 
be used in spintronics for information storage and processing. However, Skyrmions in ferromagnets 
have some intrinsic difficulties which must be overcome to use them for spintronic applications, such 
as the inability to move straight along current. We show that Skyrmions can also be stabilized 
and manipulated in antiferromagnetic materials. An antiferromagnetic Skyrmion is a compound 
topological object with a similar but of opposite sign spin texture on each sublattice, which e.g. 
results in a complete cancelation of the Magnus force. We find that the composite nature of anti¬ 
ferromagnetic Skyrmions gives rise to different dynamical behavior, both due to an applied current 
and temperature effects. 


Skyrmions are topologically protected objects which 
can form in magnetic materials [1]. They are an active 
area of research in spintronics because of their potential 
for encoding, transmitting [2, 3] and computing informa¬ 
tion [4]. Belavin and Polyakov [5] introduced Skyrmions 
[6] in the context of the two-dimensional Heisenberg 
model although Feldtkeller [7] and Thiele [8] discussed 
similar topological defects earlier in a more general con¬ 
text. Skyrmions have some advantages over other pro¬ 
posed storage technologies, such as domain wall registers, 
because of the low currents required to move them due 
to the ability to move past pinning sites [9]. Implement¬ 
ing Skyrmion devices in ferromagnetic materials involves 
issues in common with other spintronic concepts, such 
as the sensitivity to stray fields. Skyrmions also posses 
a further complication in that they experience a Mag¬ 
nus force perpendicular to the applied current, making 
it difficult to move Skyrmions along the current [10]. By 
contrast, antiferromagnets are not sensitive to stray fields 
and with an applied current we find that Skyrmions in 
antiferromagnets move in straight lines along the current, 
distinctly different from ferromagnetic materials. In an 
antiferromagnet the Skyrmion forms as a pair of strongly 
coupled topological objects, one pertaining to each sub¬ 
lattice. The opposing topological index of each sublat¬ 
tice causes an exact cancellation of the Magnus force, 
hence there is no transverse component of the velocity. 
The current induced longitudinal velocity is also found 
to strongly depend on the material parameters (a, (3) 
and as a result can reach high velocities of the order of 
km/s. Moreover, the thermal properties of antiferromag¬ 
netic Skyrmions are found to be rather different from 
their ferromagnetic counterparts. 

A compelling reason to study Skyrmions in antiferro¬ 
magnets apart from the insensitivity to stray fields, is 
that the Dzyaloshinskii-Moriya interaction (DMI), which 
is essential for the formation of individual Skyrmions, is 


more commonly found in antiferromagnetic (AFM) ma¬ 
terials than ferromagnetic (FM) materials. Most recent 
experimental results on FM Skyrmions rely on the pres¬ 
ence of an interfacial DMI to stabilize Skyrmions, how¬ 
ever bulk DMI is more prevalent in AFMs [11, 12]. AFMs 
are also considerably more abundant in nature than fer¬ 
romagnets, although metallic AFMs are not so common 
but examples include FePt 3 and M^Au. 

Skyrmions can form in different systems where there 
is a competition between the DMI and another energy 
contribution, for example the Zeeman energy from an 
applied field or a uniaxial anisotropy [3, 13]. Here the 
last option is studied by necessity because of the antifer¬ 
romagnets insensitivity to applied fields and the lack of 
a significant demagnetizing field precludes these mecha¬ 
nisms from forming a Skyrmion [14]. We also focus on 
individual Skyrmions, rather than a Skyrmion lattice, as 
the ability to move and manipulate individual bits of in¬ 
formation is more relevant to the suggested technological 
applications [2]. 

In this Letter we consider the so-called 4 G-type’ anti¬ 
ferromagnet, formed by a three-dimensional chess board 
like pattern. The AFM Skyrmion forms in much the 
same way as a FM Skyrmion, by introducing a topolog¬ 
ical defect, reversing the A and B sublattices within a 
small area and allowing the system to relax. The DMI 
prevents the metastable domain from reversing. The spin 
structure shown in Fig. 1 is analogous to the ‘hedgehog’ 
Skyrmion state of a FM but with one of the sublattices 
inverted. Hence, the topological defect exists in the Neel 
field and the magnetization is nearly zero everywhere. At 
the center of the AFM Skyrmion, neither sublattice dom¬ 
inates, but instead there is a compensation of opposite 
spins around the true center of the Skyrmion. We find 
that in the absence of temperature the radial profile and 
Skyrmion radius for a given DMI are the same in both 
FM and AFM Skyrmions for the magnetization and Neel 
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FIG. 1. The spin texture of G-type AFM Skyrmion. a) top 
view of the Skyrmion, white lines show contours of constant 
n z . The radius is 2.1 nm. b) cross-section of the Skyrmion. 
The core is not a single spin but a compensated structure 
combining the two sublattices. 


parameters respectively (Fig. 4). 

We first study the athermal dynamics of the AFM 
Skyrmion with an applied current, comparing the funda¬ 
mentals of AFM Skyrmion dynamics with Thiele’s equa¬ 
tions [8, 9, 15-17] for FM spin textures, before moving 
onto more complicated effects introduced by tempera¬ 
ture. Coupling to the current assumes that the elec¬ 
trons of up and down spin are transported predomi¬ 
nantly through their corresponding magnetization sub¬ 
lattice [18]. For the G-type AFM this is reasonable, but 
for other AFMs the transport of the electronic current 
through the AFM may be different. From this assump¬ 
tion the spatial derivative VM is calculated for the mag¬ 
netization of each sublattice, rather than the net local 
magnetization which is almost zero. 

Comparing the AFM dynamics with those of a FM 
Skyrmion (where the only change in material parameters 
is in the sign of the exchange interaction), highlights two 
main intrinsic differences in the dynamics resulting from 
the AFM characteristics. Firstly, the AFM Skyrmion al¬ 
ways has zero transverse velocity v±, relative to the cur¬ 
rent. In the FM this is only true for the highly symmetric 
case of a = /?, where a is the Gilbert damping constant 
and (3 is the non-adiabatic spin-transfer torque parame¬ 
ter. In a FM the transverse velocity is due to the Magnus 
force acting on the Skyrmion and the direction (±y) is 


determined by the winding number of the Skyrmion: 

Q {k) = J tij e Q /3 7 rn£ fc) dim^ djtn ^, (1) 

where (r) is the unit vector parallel to the local mag¬ 
netisation M^(r) and k = 1,2 label the sublattices in 
the AFM case. The AFM Skyrmion is essentially com¬ 
posed of two topological objects with opposite winding 
numbers ( Q^ = ±1) which are strongly coupled through 
the AFM exchange interaction. Both sublattices gener¬ 
ate a Magnus force, but there is a perfect cancellation 
(Fig. 2c) thus resulting in no v±. As a result the AFM 
Skyrmion travels in a perfectly straight trajectory along 
the current (see Supplementary Movie SI). One can also 
directly define the winding number for the AFM order 
parameter (Neel field) n(r, t) = m^^r, t) — m (2) M) 
and thus show that AFM Skyrmions are topologically 
nontrivial textures with the AFM topological charge ±1. 

The second notable difference between AFM and FM 
Skyrmions is that the longitudinal velocity in the AFM 
can greatly exceed the FM Skyrmion drift velocity which 
is always close to the electron drift velocity (v ~ 200 m/s 
for current j = 200 m/s and <a, /3 1). For low a or high 

/?, AFM Skyrmions can move at km/s whilst remaining 
stable (Fig. 2a). Recent theoretical studies of AFM dy¬ 
namics give an insight into this [19-21]. The dynamics 
can be studied based on the generalized Thiele’s equa¬ 
tions [19] which describe the motion of a spin texture in 
terms of collective coordinates bj 

M ij bj + aT^bj = F\ (2) 

which correspond to the soft modes of the Skyrmion [27]. 
Here AF- 7 is the mass tensor, aT^ characterizes viscous 
friction and is related to damping in the AFM, and fi¬ 
nally F l is the generalized force due to the current. The 
mass term contributes to the dynamics only on short 
time scales (up to ^ 2 ps), as shown in the inset of 
Fig. 2a. The AFM Skyrmion reaches its terminal ve¬ 
locity quickly, which is understood more intuitively by 
switching back to the two-sublattice description and writ¬ 
ing the Thiele’s equations for m( 1,2 )( r , t) [27]. It is then 
clear that the Magnus forces for (r, t) and m (2) M), 
Gzxj with the gyrocoupling constant G = 47 xQ^ k \ cancel 
each other (Fig. 2c) and the remaining generalized drag 
force, F = T/3j, leads to the AFM Skyrmion velocity 

P • (o\ 

^|| = -J 3 

11 a 

which is only along the direction of the current. The 
velocity is plotted as lines in Fig. 2a, showing an excellent 
agreement with the simulations. 

The importance of understanding the thermal prop¬ 
erties of Skyrmions is now becoming clear [22]. On a 
macroscopic level Skyrmions diffuse due to the thermal 
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FIG. 3. a) Brownian motion of the AFM Skyrmion for Gilbert 
damping constants a = 0.1 and 0.01 at T/T c — 0.25. b) 
Diffusion coefficient of the AFM Skyrmion as a function of 
a. Points are calculated from the mean squared displacement 
(r 2 ) (inset), the solid line is V — XksT / (2asa). 


FIG. 2. Current induced AFM Skyrmion dynamics (j — 200 
m/s), a) longitudinal velocity for different combinations of 
a and /3. Points are calculated numerically and the lines are 
Eq. (3) based on Thiele’s equations for an AFM. Inset shows 
the mass term is small and the Skyrmion reaches terminal ve¬ 
locity after 2 ps. b) transverse velocities calculated from the 
same simulations show there is no transverse motion, c) the 
AFM Skyrmion is composed of two topological objects with 
opposite topological charge, hence the Magnus force acts in 
opposite directions. The strong coupling between the sublat¬ 
tices leads to a perfect cancellation of the two opposing forces 
and so the AFM Skyrmion has no transverse motion. 

perturbations of the magnetic moments. We also con¬ 
sider the temperature dependence of the macroscopic ma¬ 
terial parameters, such as anisotropy and exchange stiff¬ 
ness, and find this can lead to a change in the balance 
of the competing energy terms in Skyrmionic systems. 
The Skyrmion spin texture is also subjected to deforma¬ 
tions due to internal dynamics which are stimulated by 
thermally induced spin waves. Using Langevin Landau- 
Lifshitz-Gilbert simulations [27] we have studied these 
thermal effects for the AFM Skyrmion and compared 
it with that of the FM Skyrmion (see Supplementary 


Movies S3 and S4). 

Simulating the Brownian motion of a single Skyrmion 
(Fig. 3), we find the AFM Skyrmion to be diffusive, 
meaning the mean square displacement (r 2 ) oc £, as was 
shown for the FM Skyrmion [22]. However the diffusion 
coefficient of the AFM Skyrmion is greater than that of 
the FM and D oc 1/ct as a direct consequence of the 
absence of a net Magnus force which plays the dominant 
role in the FM diffusion where D oc a [22]. In Ref. 23 the 
diffusion coefficient of AFM textures in one-dimensional 
case was shown to be V = XksT/(2ascr), i.e. inversely 
proportional to the Gilbert damping <a, where s = hS/a 3 , 
A = A/K is the exchange length and a is the cross 
sectional area of the domain wall - in this case the cir¬ 
cumferential area of the Skyrmion, 2iraR s (T). The lack 
of a Magnus force allows this expression to be general¬ 
ized to two-dimensional diffusion, we compared it with 
our numerical results (solid line in Fig. 3, no fitting was 
performed) finding a good agreement. 

The increased thermal mobility of the AFM Skyrmions 
may be useful if attempting to move spin textures with 
heat gradients. However, it also poses challenges for more 
conventional current driven motion as the thermal per- 


































4 



FIG. 4. AFM and FM Skyrmion radii R s as a function of 
reduced temperature T/Tc- The error bars represent the 
magnitude of the thermally induced radius fluctuations. The 
solid line is the simple scaling theory, Eq. (4), dashed lines 
are power law fits R S (T) oc ra -0 ' 86 and R S (T) oc n -2 ' 09 . 


turbations may cause excessive randomness in the mo¬ 
tion. We have studied just one of the AFM types in 
this Letter and it is thus possible that other AFMs could 
contain the benefit of high drift velocity but with a lower 
diffusion coefficient. 

We also calculated the temperature dependence of the 
mean Skyrmion radius, R s , and compared to an estimate 
from the scaling relationships of A , D and K with the 
reduced magnetization m = 0). The scaling of 

A oc m 2 and K oc ra 3 in bulk FMs are well known and 
we use these as approximations in the thin film. In the 
absence of knowledge concerning the thermal scaling of 
D we find that assuming no temperature dependence of 
D gives the closest agreement with the numerical results 
(Fig. 4), the final equation being 


Rs(T) = 


2AX 


4y/AKm(T) - irDm-^/ 2 {T )' 


( 4 ) 


The results show that increasing temperature causes 
a change in the domain wall energy cost, with the ef¬ 
fective DMI contribution becoming larger as tempera¬ 
ture increases, causing an increase in the Skyrmion ra¬ 
dius [24]. Surprisingly the AFM Skyrmion radius shows 
a stronger temperature dependence and larger fluctua¬ 
tions of the radius about the mean value than for the 
FM Skyrmion. At zero temperature we found the FM 
and AFM Skyrmion to have the same R s vs. D depen¬ 
dence as shown by Rohart et al. [25] for FM Skyrmions. 
Hence the deviation observed at finite temperature is the 
result of the different characteristic thermal fluctuations 
in an FM and AFM. 

The plethora of different types of antiferromagnetic 
materials makes further research in the area of AFM 
Skyrmions an exciting prospect. Searching for AFM 
Skyrmions experimentally will be a challenging task, but 


techniques such as neutron scattering have the potential 
to find AFM Skyrmion lattices. Moreover, x-ray mag¬ 
netic linear dichroism occurs in the AFM state as the 
spin-orbit coupling leads to a distortion of the charge 
density and, thus the AFM structure can be experimen¬ 
tally measured [26]. 

Small spin textures can serve as bits of information and 
manipulating them by electric or thermal currents one of 
the main challenges in the field of spintronics. Ferro¬ 
magnetic Skyrmions recently attracted a lot of attention 
because of their small size and ability to avoid pinning 
while moved by electric current better than domain walls. 
However, ferromagnetic Skyrmions still suffer from the 
detrimental effects of stray fields and transverse intrinsic 
dynamics causing difficulties in employing them in spin- 
tronic applications. The related topological object we 
explored here - the AFM Skyrmion - overcomes these 
disadvantages, having no demagnetizing field, an insensi¬ 
tivity to stray fields, and we have shown that its dynam¬ 
ics are strictly along the current (no Hall effect) while 
potentially being faster compared to its ferromagnetic 
analogue. This makes AFM Skyrmions an ideal informa¬ 
tion carrier. However, the thermal properties, such as the 
AFM Skyrmion radius and diffusion constant differ from 
those for ferromagnetic Skyrmions because of the funda¬ 
mental difference in the spin correlations. More studies 
into the different AFM types and the effect of the spin 
correlations are needed to discover if the benefits of the 
AFM Skyrmions can be achieved but with a reduced dif- 
fusivity than we have found here. 
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Note added in proof. - Recently, we became aware of a 
recent report [30] on another study of AFM Skyrmions. 
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